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ABSTRACT 

An important concern in the application of gamma-ray bursts (GRBs) to cosmology is that the calibration 
of GRB luminosity/energy relations depends on the cosmological model, due to the lack of a sufficient low- 
redshift GRB sample. In this paper, we present a new method to calibrate GRB relations in a cosmology- 
independent way. Since objects at the same redshift should have the same luminosity distance and since the 
distance moduli of Type la supernovae (SNe la) obtained directly from observations are completely cosmology 
independent, we obtain the distance modulus of a GRB at a given redshift by interpolating from the Hubble 
diagram of SNe la. Then we calibrate seven GRB relations without assuming a particular cosmological model 
and construct a GRB Hubble diagram to constrain cosmological parameters. From the 42 GRBs at 1.4 < z < 
6.6, we obtain M = 0.25^;^, = 0.75+JJ H for the flat ACDM model, and for the dark energy model with a 
constant equation of state wq = — 1.05^^, which is consistent with the concordance model in a l-a confidence 
region. 

Subject headings: gamma rays: bursts — cosmology: observations 



1. INTRODUCTION 

In the past decade, observations of Type la supernovae (SNe 
la; Riess et al. 1998; Schmidt et al. 1998; Perlmutter et 
al. 1999), cosmic microwave background (CMB) fluctuations 
(Bennett et al. 2003; Spergel et al. 2003, 2007), and large- 
scale structures (LSS; Tegmark et al. 2004, 2006) have been 
used to explore cosmology extensively. These observations 
are found to be consistent with the so-called concordance cos- 
mology, in which the universe is spatially flat and contains 
pressureless matter and dark energy, with fractional energy 
densities of fl M = 0.27 ± 0.04 and Q, A = 0.73 ± 0.04 (Davis et 
al. 2007). Observations of SN la provide a powerful probe in 
the modern cosmology. Phillips (1993) found that there is an 
intrinsic relation between the peak luminosity and the shape 
of the light curve of SNe la. This relation and other similar 
luminosity relations make SNe la standard candles for mea- 
suring the geometry and dynamics of the universe. However, 
the maximum redshift of the SNe la which we can currently 
use is only about 1.7, whereas fluctuations of the CMB pro- 
vide the cosmological information from last scattering surface 
at z = 1089. Therefore, the earlier universe at higher redshift 
may not be well-studied without data from standardized can- 
dles in the "cosmological desert" from the SNe la redshift 
limit to z~ 1000. 

Gamma-Ray Bursts (GRBs) are the most intense explosions 
observed so far. Their high energy photons in the gamma- 
ray band are almost immune to dust extinction, whereas in 
the case of SN la observations, there is extinction from the 
interstellar medium when optical photons propagate towards 
us. Moreover, GRBs are likely to occur in the high-redshift 
range up to at least z = 6.6 (Rrimm et al. 2006); higher red- 
shift GRBs up to z = 10 should have already been detected, 
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although none have been identified (Lamb & Reichart 2000; 
Bromm & Loeb 2002, 2006; Lin, Zhang, & Li 2004). Thus, 
by using GRBs, we may explore the early universe in the high 
redshift range which is difficult to access by other cosmolog- 
ical probes. These advantages make GRBs attractive for cos- 
mology research. 

GRB luminosity/energy relations are connections between 
measurable properties of the prompt gamma-ray emission and 
the luminosity or energy. In recent years, several empirical 
GRB luminosity relations have been proposed as distance in- 
dicators (see e.g. Ghirlanda et al. 2006a; Schaefer (2007) for 
reviews), such as the luminosity-spectral lag (L-Ti ag ) relation 
(Norris, Marani, & Bonnell 2000), the luminosity-variability 
(L-V) relation (Fenimore & Ramirez-Ruiz 2000; Reichart et 
al. 2001), the isotropic energy-peak spectral energy ((EW,- 
E p ) relation (i.e., the so-called Amati relation, Amati et al. 
2002), the collimation-corrected energy-peak spectral energy 
(E^-Ep) relation (i.e., the so-called Ghirlanda et al. 2004a), 
the L-£p relation (Schaefer 2003a; Yonetoku et al. 2004), and 
two multi-variable relations. The first of these multiple re- 
lations is between E[ so , E p and the break time of the optical 
afterglow light curves (fb) (i.e., the so-called Liang-Zhang re- 
lation, Liang & Zhang 2005); the other is between luminosity, 
E p and the rest-frame "high-signal" timescale (ro.45) (Firmani 
et al. 2006). 

Several authors have made use of these GRB luminosity 
indicators as standard candles at very high redshift for cos- 
mology research (e.g. Schaefer 2003b; Dai et al. 2004 
Ghirlanda et al. 2004b; Firmani et al. 2005, 2006, 2007 
Liang & Zhang 2005; Xu et al. 2005; Wang & Dai 2006 
see also e.g. Ghirlanda et al. 2006a and Schaefer (2007) for 
reviews). Schaefer (2003b) derived the luminosity distances 
of nine GRBs with known redshifts by using two GRB lumi- 
nosity relations to construct the first GRB Hubble diagram. 
Dai et al. (2004) considered the Ghirlanda relation with 12 
bursts and proposed another approach to constrain cosmolog- 
ical parameters. Liang & Zhang (2005) constrained cosmo- 
logical parameters and the transition redshift using the E lso - 
Ep-t\, relation. More recently, Schaefer (2007) used five GRB 
relations calibrated with 69 GRBs by assuming two adopted 
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cosmological models to obtain the derived distance moduli 
for plotting the Hubble diagram, and joint constraints on the 
cosmological parameters and dark energy models have been 
derived in many works by combining the 69 GRB data with 
SNe la and the other cosmological probes, such as the CMB 
anisotropy, the baryon acoustic oscillation (BAO) peak, the X- 
ray gas mass fraction in clusters, the linear growth rate of per- 
turbations, and the angular diameter distances to radio galax- 
ies (Wright 2007; Wang et al. 2007; Li et al. 2008; Qi et al. 
2008; Daly et al. 2008). 

However, an important point related to the use of GRBs for 
cosmology is the dependence on the cosmological model in 
the calibration of GRB relations. In the case of SN la cos- 
mology, the calibration is carried out with a sample of SNe 
la at very low redshift where the luminosities of SNe la are 
essentially independent of any cosmological model (i.e., at 
Z < 0.1, the luminosity distance has a negligible dependence 
on the choice of the model). However, in the case of GRBs, 
the observed long-GRB rate falls off rapidly at low redshifts, 
and some nearby GRBs may be intrinsically different (e.g., 
GRB 980425, GRB 031203; Norris 2002; Soderberg et al. 
2004; Guetta et al. 2004; Liang & Zhang 2006a). Therefore, 
it is very difficult to calibrate the relations with a low-redshift 
sample. The relations of GRBs presented above have been 
calibrated by assuming a particular cosmological model (e.g. 
the ACDM model). In order to investigate cosmology, the 
relations of standard candles should be calibrated in a cos- 
mological model-independent way. Otherwise, the circularity 
problem cannot easily be avoided. Many of the works men- 
tioned above treat the circularity problem with a statistical 
approach. A simultaneous fit of the parameters in the cali- 
bration curves and the cosmology is carried out to find the 
optimal GRB relation and the optimal cosmological model 
in the sense of a minimum scattering in both the luminos- 
ity relations and the Hubble diagram. Firmani et al. (2005) 
has also proposed a Bayesian method to get around the cir- 
cularity problem. Li et al. (2008b) presented another new 
method to deal with the problem, using Markov chain Monte 
Carlo (MCMC) global fitting analysis. Instead of a hybrid 
sample over the whole redshift range of GRBs, Takahashi et 
al. (2003) first calibrated two GRB relations at low redshift 
(z < 1), where distance-redshift relations have been already 
determined from SN la; they then used GRBs at high redshift 
(z > 1) as a distance indicator. Bertolami & Silva (2006) con- 
sidered the use of GRBs at 1 .5 < z < 5 calibrated with the 
bursts at z < 1 .5 as distance markers to study the unification 
of dark energy and dark matter in the context of the general- 
ized Chaplygin gas model. Schaefer (2007) made a detailed 
comparison between the Hubble diagram calibrated only with 
37 bursts at z < 2 and the one calibrated with the data from all 
69 GRBs to show that the results from fits to the GRB Hubble 
diagram calibrated with the hybrid sample are robust. These 
calibration methods above carried out with the sample at low 
redshifts have been derived from the ACDM model. However, 
we note that the circularity problem can not be circumvented 
completely by means of statistical approaches, because a par- 
ticular cosmology model is required in doing the joint fitting. 
This means that the parameters of the calibrated relations are 
still coupled to the cosmological parameters derived from a 
given cosmological model. In principle, the circularity prob- 
lem can be avoided in two ways (Ghirlanda et al. 2006a): 
(1) through a solid physical interpretation of these relations 
that would fix their slope independently from cosmological 
model, or (2) through the calibration of these relations by sev- 



eral low-redshift GRBs. Recently the possibility of calibrat- 
ing the standard candles using GRBs in a low-dispersion in 
redshift near a fiducial redshift has been proposed (Lamb et al. 
2005; Ghirlanda et al. 2006b). Liang & Zhang (2006b) elabo- 
rated this method further based on Bayesian theory. However, 
the GRB sample available now is far from what is needed to 
calibrate the relations in this way. 

As analyzed above, due to the lack of low-redshift GRBs, 
the methods used in GRB luminosity relations for cosmol- 
ogy are different from the standard Hubble diagram method 
used in SN la cosmology. In this work, we present a new 
method to calibrate the GRB relations in a cosmological 
model-independent way. In the case of SN la cosmology, 
the distance of nearby SNe la used to calibrate the luminos- 
ity relations can be obtained by measuring Cepheid variables 
in the same galaxy, and with other distance indicators. Thus, 
Cepheid variables have been regarded as the first-order stan- 
dard candles for calibrating SNe la, the second-order standard 
candles. It is obvious that objects at the same redshift should 
have the same luminosity distance in any cosmology. For cal- 
ibrating GRBs, there are so many SNe la (e.g., 192 SNe la 
used in Davis et al. 2007) that we can obtain the distance 
moduli (also the luminosity distance) at any redshift in the 
redshift range of SNe la by interpolating from SN la data in 
the Hubble diagram. Furthermore, the distance moduli of SNe 
la obtained directly from observations are completely cosmo- 
logical model independent. Therefore, in the same sense as 
Cepheid variables and SNe la, if we regard SNe la as the first 
order standard candles, we can obtain the distance moduli of 
GRBs in the redshift range of SNe la and calibrate GRB rela- 
tions in a completely cosmological model-independent way. 
Then if we further assume that these calibrated GRB relations 
are still valid in the whole redshift range of GRBs, just as for 
SNe la, we can use the standard Hubble diagram method to 
constrain the cosmological parameters from the GRB data at 
high redshift obtained by utilizing the relations. 

The structure of this paper is arranged as follow. In section 
2, we calibrate seven GRB luminosity/energy relations with 
the sample at z < 1 .4 obtained by interpolating from SNe la 
data in the Hubble diagram. In section 3, we construct the 
Hubble diagram of GRBs obtained by using the interpolation 
methods and constrain cosmological parameters. Conclusions 
and a discussion are given in section 4. 

2. THE CALIBRATION OF THE LUMINOSITY RELATIONS OF 
GAMMA-RAY BURSTS 

We adopt the data for 192 SNe la (Riess et al. 2007; Wood- 
Vasey et al. 2007; Astier et al. 2007, Davis et al. 2007), 
shown in Figure 1. It is clear that for GRBs in the redshift 
range of SNe la, there are enough data SN la points and the 
redshift intervals of the neighboring SN la data points are also 
small enough, to be able to use interpolation methods to ob- 
tain the distance moduli of GRBs. Since there is only one 
SN la point at z > 1.4 (the redshift of SN1997ff is z = 1.755), 
we initially exclude it from our SN la sample used to inter- 
polate the distance moduli of GRBs in the redshift range of 
the SN la sample (191 SNe la data at z < 1.4). As the dis- 
tance moduli of SNe la data in the Hubble diagram are ob- 
tained directly from observations, with the sample at z < 1.4 
by interpolating from the Hubble diagram of SNe la, we can 
therefore calibrate GRB luminosity relations in a completely 
cosmology-independent way. 

In this section, we calibrate seven GRB luminosity/energy 
relations with the sample at z < 1 .4, i.e., the ri ag -L relation, the 
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V-L relation, the L-E p relation, the E 1 -E v relation, the trt-L 
relation, where trt is the minimum rise time in the GRB light 
curve (Schaefer, 2007), the E lso -E p relation, and the E-^o-Ep-tb 
relation. 

The isotropic luminosity of a burst is calculated by 

L = 4irdlP boio , (1) 

where d L is the luminosity distance of the burst and fb ] is 
the bolometric flux of gamma-rays in the burst. The isotropic 
energy released from a burst is given by 

E h0 =47r^S bo io(l +zr\ (2) 

where 5boio lS the bolometric fluence of gamma-rays in the 
burst at redshift z. The total collimation-corrected energy is 
then calculated by 

— ^beam^iso: (3) 

where the beaming factor, Fb eam is (l-cos#j et ) with the jet 
opening angle (0j e t), which is related to the break time (?b). 

A GRB luminosity relation can be generally written in the 
form 

logy = a + b\ogx, (4) 

where a and b are the intercept and slope of the relation re- 
spectively; y is the luminosity (L in units of erg s" 1 ) or energy 
(E i so or E 1 in units of erg); x is the GRB parameters mea- 
sured in the rest frame, e.g., T lag (l+z) _1 /(0.1 s), V(l+z)/0.02, 
E p (l +z)/(300 keV), r RT ( 1 +zf 1 /(0. 1 s), £ p (l +z)/(300 keV), 
for the 6 two-variable relations above. We adopt the data 
for these quantities from Schaefer (2007). For the one multi- 
variable relation (i.e., Ei so -E p -t\,), the calibration equation is 

logy = a + Mogxi+b 2 log.*2, (5) 

where x\ and x 2 are £ p (l +z)/(300 keV), f b /(l +z)/(l day) 
respectively, and b\ and b 2 are the slopes of x\ and x 2 respec- 
tively. The error of the interpolated distance modulus of a 
GRB must account for the original uncertainties of the SNe la 
as well as for the uncertainties from the interpolation. When 
the linear interpolation is used, the error can be calculated by 

o» = ([(Zi+i -z)/(z;+i -zdfe 2 i+[(z-zd/(zm -zdfe 2 ,- +1 ) 1/2 , 

(6) 

where is the error of the interpolated distance modulus, 
/i is the interpolated distance modulus of a source at redshift 
Z, €fj,j and e M ,/+i are errors of the SNe, /x,- and tt,+i are the 
distance moduli of the SNe at nearby redshifts Zi and Zi+i re- 
spectively. When the cubic interpolation is used, the error can 
be calculated by 

°> = (^0 e ^,/ + ^l e ^,!+l + ^2 e ^,!+2 + ^3 e /V+3) 1//2 i (J) 

where e^j, e ft ;+i, e^j+ 2 , and e^j+3 are errors of the SNe; and 
(ttj, fii+i, iij+2, and /x,+3 are the distance moduli of the SNe at 
nearby redshifts zu Zi+i, Zi+ 2 , and z i+3 : 

Ao = [(Zi+l -Z)(Z»+2-Z)fe+3-Z)]/[(Zi+l -Zi)(Zi+ 2 -Zi)(Zi+3-Zi)V, 
Ai = [(Zi~ Z)(Z/+2-Z)(Zi+3-Z)]/[(Z/- Zi+l)(Z/+2-Zi+l)(Z/+3-Zi+l)]; 
M = [(Zi-Z)(Z !+ l -Z)(Zi+3-Z)]/[(Zi- Zi+2)(Zi+l -Zi+2)(Zi+3-Zi+2)]; 
^3 = [(Z,— Z)(Zi+l -Z)(Zi+2-Z)]/[(Zi-Zi+3)(Z(+l -Zm){Zi+2-Zi+l)~\. 

We determine the values of the intercept (a) and the slope 
(b) with their 1-er uncertainties calibrated with the GRB sam- 
ple at z < 1 .4 by using two interpolation methods (the linear 
interpolation methods and the cubic interpolation method). 
For the 6 two-variable relations we use the same method (the 
bisector of the two ordinary least-squares) as used in Schae- 
fer (2007), and for the one multi-variable relation, the multi- 
ple variable regression analysis is used. The bisector of the 
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FIG. 1 .— Hubble Diagram of 192 SNe la (red dots) and the 69 GBRs (cir- 
cles) obtained using the interpolation methods. The 27 GBRs at z < 1.4 are 
obtained by interpolating from SN la data (black circles: the cubic interpola- 
tion method; black stars: the linear interpolation method); and the 42 GBRs 
at z > 1 .4 (blue circles) are obtained with the five relations calibrated with 
the sample at z < 1.4 using the cubic interpolation method. The curve is the 
theoretical distance modulus in the concordance model (w>o = —1, Qm = 0.27, 
Ca = 0.73), and the vertical dotted line represents z = 1.4. 

two ordinary least-squares (Isobe et al., 1990) does not take 
the errors into account; but the use of weighted least-squares, 
taking into account the measurable uncertainties, results in al- 
most identical best fits. However, taking into account the mea- 
surable uncertainties in the regression, when they are smaller 
than the intrinsic error, indeed will not change the fitting pa- 
rameters significantly (for further discussion, see Schaefer 
2007). The calibration results are summarized in Table 1, and 
we plot the GRB data at z < 1.4 with the distance moduli 
obtained by using the two interpolation methods from SN la 
data in Figure 1 for comparison. In the previous treatment, 
the distance moduli of GRBs are obtained by assuming a par- 
ticular cosmological model with a hybrid sample in the whole 
redshift range of GRBs to calibrate the relations. Therefore, 
for comparison in Table 1 we also list the results calibrated 
with the same sample (z < 1 .4) assuming the ACDM model 
(w = -1) or the Riess cosmology 4 (w = —1.31 + 1.48z, Riess et 
al. 2004), which were used to calibrate the five GRB relations 
in Schaefer (2007). The calibration results obtained by using 
the interpolation methods are carried out with the 27 GRBs at 
z < 1.4, namely, 13, 19, 25, 12, 24, 12, and 12 GRBs for the 
Ti ag -L, V-L, L-E p , E-f-Ep, trt-L, ^-^p, and E iti0 -E p -t b rela- 
tions respectively; 5 therefore, the uncertainties of the results 
are somewhat larger than those with the hybrid GRB sample. 6 



Here we follow Schaefer (2007) in calling this particular parametrization 
of the equation of state the "Riess cosmology". 

5 For the seven relations, only data from Table 4 in Schaefer (2007) are 
used for the calibration. However, we note that there are more GRB data 
available to calibrate the Ei s0 -Ep relation. 

6 Calibration results with their 1-cr uncertainty for the five GRB relations 
with 69 GRBs (0. 1 < z < 6.6) obtained by assuming the ACDM model in 
Schaefer (2007) are also given in Table 1 , for comparison. 
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table l 
Calibration results 



Relation 


Linear interpolation method 


Cubic interpolation method 


The ACDM model 


The Riess 


cosmology 


Schaefer (2007) 


a 


b 


a 


b 


a 


b 


a 


b 


a 


b 


T)ag-£ 


52.22±0.09 


-1.07±0.14 


52.22±0.09 


-1.07±0.13 


52. 15 ±0.10 


-1.11±0.14 


52. 13 ±0.09 


-1.10±0.13 


52.26±0.06 


-1.01±0.05 


V-L 


52.58±0.13 


2.04±0.26 


52.59±0.13 


2.05±0.27 


52.48±0.13 


2.04±0.30 


52.47±0.13 


2.02±0.30 


52.49±0.22 


1.77±0.12 


L-E ¥ 


52.25 ±0.09 


1.68±0.11 


52.26±0.09 


1.69±0.11 


52. 15 ±0.09 


1.66±0.12 


52. 15 ±0.09 


1.64±0.12 


52.21 ±0.13 


1.68±0.05 




50.70±0.06 


1.77±0.17 


50.71 ±0.07 


1.79±0.18 


50.59±0.06 


1.75±0.19 


50.59±0.06 


1.73±0.19 


50.57±0.09 


1.63 ±0.03 


T RT -L 


52.63 ±0.11 


-1.31±0.15 


52.64±0.11 


-1.31±0.15 


52.52±0.10 


-1.30±0.15 


52.51 ±0.10 


-1.29±0.15 


52.54±0.06 


-1.21 ±0.06 




52.99±0.16 


1.88±0.23 


53.00±0.17 


1.90±0.23 


52.88±0.16 


1.84±0.21 


52.88±0.15 


1.81±0.20 








52.83±0.10 


2.26±0.30 


52.83±0.10 


2.28±0.30 


52.73 ±0.10 


2.19±0.30 


52.73 ±0.10 


2.16±0.30 










-1.07±0.21 




-1.07±0.21 




-1.01±0.20 




-0.99±0.21 







NOTE. — Calibration results (for o=intercept, £>=slope) with their 1-cr uncertainties, for the seven GRB luminosity/energy relations with the 
sample at z < 1 .4, using two interpolation methods (the linear interpolation method and the cubic interpolation methods) directly from SNe la 
data and by assuming a particular cosmological models (the ACDM model or the Riess cosmology). Only the data available from Table 4 in 
Schaefer 2007 are used to calibrate the seven GRB relations. For comparison, calibration results with the 69 GRBs obtained by assuming the 
ACDM model for the five GRB relations in Schaefer 2007 are given in the last column. 



"In the slope values for the Ei so -E p -tb relation, the first line is b\, and the 
second line is b2- 
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The linear correlation coefficients of the calibration with 
the sample at z < 1.4 using the cubic interpolation methods 
are -0.88, 0.65, 0.89, 0.94, -0.75, and 0.68 for the above 6 
two-variable relations respectively, and 0.94 for E iso vs. the 
combined variable (b\ log[2? p (l+z)]+&2log|/b/(l +z)]) f° r the 
£iso-£p-fb relation, which shows that these correlations are 
significant. 

From Table 1 and Figure 1, we find that the calibration re- 
sults obtained using the linear interpolation methods are al- 
most identical to the results calibrated by using the cubic in- 
terpolation method. We also find the results obtained by as- 
suming the two cosmological models with the same sample 
differ only slightly from, but are still fully consistent with, 
those calibrated using our interpolation methods. The reason 
for this is easy to understand, since both cosmological models 
are fully compatible with SN la data. Nevertheless, it should 
be noticed that the calibration results obtained using the in- 
terpolation methods directly from SN la data are completely 
cosmology independent. 

3. THE HUBBLE DIAGRAM OF GAMMA-RAY BURSTS 

If we further assum GRB luminosity/energy relations do not 
evolve with redshift, we are able to obtain the luminosity (L) 
or energy (£i S0 or E 7 ) of each burst at high redshift (z > 1 .4) 
by utilizing the calibrated relations. Therefore, the luminosity 
distance (<f L ) can be derived from equation (1~3). The uncer- 
tainty of the value of the luminosity or energy deduced from 
a GRB relation is 



cLy = °l + (<rb\ogx) 2 + (0.4343ba x /x) 2 + a l 



2 

sys ' 



(8) 



where a a , Ob and a x are l-cr uncertainty of the intercept, the 
slope and the GRB measurable parameters, and er sys is the sys- 
tematic error in the fitting that accounts for the extra scatter of 
the luminosity relations. The value of er sys can be estimated 
by finding the value such that a x 2 fit to the calibration curve 
produces a value of reduced x 2 of unity (Schaefer 2007). A 
distance modulus can be calculated as fi = 51ogc/ L +25 (where 
d L is in Mpc). The propagated uncertainties will depend on 
whether Pboio or Sboio is used: 



<r M = [(2.5a, ogL ) 2 + (1.086a ftolo / J P bolo ) z ] 



2-1 1/2 



(9) 



or 



and 



tr„ = [(2.5a log£iso ) 2 + (1.086 ( r Sbulo /5boio) 2 ] 1/2 , (10) 



a M = [(2.5 ( x log ^) 2 +(l .086aj bolo /5 bolo ) 2 +(l .086^/Fbe^) 2 ] l ' 2 . 

(11) 

We use the same method as used in Schaefer (2007) to ob- 
tain the best estimate /i for each GRB which is the weighted 
average of all available distance moduli. The derived distance 
modulus for each GRB is 



(12) 



with its uncertainty ct m = o-" 2 )~ x l 2 , where the summations 
run from 1 to 5 over the five relations used in Schaefer (2007) 
with available data. 

We have plotted the Hubble diagram of the 69 GRBs ob- 
tained using the interpolation methods in Figure 1. The 27 
GRBs at z < 1 4 are obtained by using two interpolation meth- 
ods directly from SNe data. The 42 GRB data at z > 1 4 are 
obtained by utilizing the five relations calibrated with the sam- 
ple at z < 1 -4 using the cubic interpolation method. 



Then the cosmological parameters can be fitted by the min- 
imum x 2 method. The definition of x 2 is 



X 



(Mtt 



"Mobs,;) 



(13) 



in which , is the theoretical value of distance modulus, 
^ b s , i is the observed distance modulus with its error . The 
theoretical value of the distance modulus (^th) depends on 
the theoretical value of luminosity distance. For the ACDM 
model, the luminosity distance is calculated by 



c(l+z) 



x sinn[-\ 



if 

Jo 



dz 



V(i+z) 2 (i+n m z)-z(2+z)n A 

(14) 

where ilk = 1 - _ ^A, and sinn(x) is sinh for ilk > 0, sin 
for ilk < 0, and x for ilk = 0. For the dark energy model with 
a constant equation of state (wq), the luminosity distance in a 
flat universe is 



, c(l+z) 
d L = — ~ — x 



dz 



Ho JO y/l 1 + z) 3 il m + ( 1 - OmX 1 + Z) 3(1+W0) 

(15) 

here we adopt Ho = 70 km s 'Mpc 1 . 

Figure 2a shows the joint confidence regions for (i! M ,ilA) 
in the ACDM model from the 42 GRB data (z > 1.4) ob- 
tained by utilizing the five relations calibrated with the sam- 
ple at z < 1-4 using the cubic interpolation method. The 1- 
a confidence region is (O m ,^a) = (0.25^,0.68;^) with 
Xmin = 44.40 for 40 degrees of freedom. For a flat universe, 
we obtain il M = 0.25^$ and il A = 0.75^$, which is con- 
sistent with the concordance model (51m = 0.27, ilA = 0.73) 
in the l-cr confidence region. From equation (14), it is ob- 
vious that for the ACDM model, the theoretical value of the 
luminosity distance mainly depends on il M at higher redshift, 
and strongly depends on ilA at lower redshift. Therefore, we 
can find that the shape of the likelihood contour of GRBs at 
higher redshift is almost vertical to the horizontal axis (iljvi) 
compared to that of SNe la at lower redshift. To understand 
the shape of the confidence regions in the (SIm^a) plane, 
Firmani et al. (2007) explored the behavior of the luminos- 
ity distance d^ at different redshifts z in a given cosmolog- 
ical parameter space. The stripe where <f L varies by 1% at 
Z = 3, shown in the (S! M , £1a) plane, was almost vertical to the 
il M axis, which corresponds roughly to the typical redshifts 
(z ~ 3) of the GRB sample at z > 1 .4. If we add the one SN 
la point (SN1997ff at z = 1.755) at z > 1.4 into our SN la 
sample used to interpolate the distance moduli of GRBs, we 
can calibrate the GRB relations with 36 GRBs at z < 1 .76 by 
using the cubic interpolation method. We find that the fitting 
results from the 33 GRB data (z > 1.76) obtained by utiliz- 
ing the relations calibrated with the sample at z < 1 .76 using 
the cubic interpolation method is SIm = 0.28^q q7 for a flat uni- 
verse (Figure 2b), which is consistent with the result from 
the 42 GRB data (z > 1 .4) obtained using the cubic interpo- 
lation method. However, if the relations are calibrated with 
the GRB sample at z < 1 .4 by assuming a particular cosmo- 
logical model, we find that the fitting results the 42 GRB data 
(z > 1 .4) obtained assuming the ACDM model or the Riess 
cosmology are M = 0.38!j]$ (Figure 2c) or M = 0.39^ 
(Figure 2d) for a flat universe, which are systematically higher 



than il M = 0.27±g Q4 beyond a l-cr deviation. 

Figure 3 shows likelihood contours in the (il M , wq) plane in 
the dark energy model with a constant wq for a flat universe 
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FIG. 2. — (a) Joint confidence regions for (Qm,^a) in the ACDM model 
from the data for 42 GRBs (z > 1.4) obtained by utilizing the five relations 
calibrated with the sample at z < 1.4 using the cubic interpolation method. 
The plus sign indicates the best fit values. The contours correspond to 1, 
2, and 3-cr confidence regions, and the dashed line represents the flat uni- 



verse. The l-cr confidence region are (Qm,^a) = (0.25 j ! 



.().( 



1). For 



a flat universe prior, Cm = 0.25 jj [g and = 0.75+jjjjJ. (b) Joint confidence 
regions for (Hm,Ca) from the data for 33 GRBs (z > 1.76) obtained by uti- 
lizing the relations calibrated with the sample at z < 1.76 by using the cubic 
interpolation method. For a fiat universe prior, f2jy] = 0.28^q q^. (c) The joint 
confidence regions for (Cm^a) from the data for 42 GRBs (z > 1.4) ob- 
tained by utilizing the relations calibrated with the sample at z < 1 -4, assum- 
ing a ACDM model. For a flat universe prior, f2jyj = 0.38^] J]^. (d) The joint 
confidence regions for (Cm^a) from the 42 GRBs data (z > 1.4) obtained 
by utilizing the relations calibrated with the sample at z < 1 -4, assuming a 
Riess cosmology. For a flat universe prior, f2jyi = 0.39^ Qy. 




FIG. 3. — Contours of likelihood in the (Qm, wo) plane in the dark energy 
model with a constant wq for a flat universe from the 42 GRBs data (z. > 1.4) 
obtained by utilizing the five relations calibrated with the sample at z < 1.4 
using the cubic interpolation method. The plus sign indicates the best-fit 
values (Qm = 0.25, wq = —0.95). The contours correspond to 1, 2, and 3-cr 
confidence regions. For a prior of f2jyi = 0.27 (vertical line), wq 
(solid part of vertical line). 



, qc+0.40 



from the 42 GRB data (z > 1 .4) obtained by utilizing the five 
relations calibrated with the sample at z < 1.4 using the cubic 
interpolation method. The best-fit values are SIm = 0.25, wq = 
-0.95. For a prior of 17m = 0.27, we obtain wq = -1.05^ ^, 
which is consistent with the cosmological constant in a l-cr 
confidence region. 



4. SUMMARY AND DISCUSSION 

With the basic assumption that objects at the same redshift 
should have the same luminosity distance, we can obtain the 
distance modulus of a GRB at given redshift by interpolating 
from the Hubble diagram of SNe la at z < 1 .4. Since the dis- 
tance modulus of SN la is completely cosmological model in- 
dependent, the GRB luminosity relations can be calibrated in 
a completely cosmology independent way. Instead of a hybrid 
sample in the whole redshift range of GRBs used in most pre- 
vious treatment, we choose the GRB sample only in the red- 
shift range of SNe la to calibrate the relations. We find there 
are not significant differences between the calibration results 
obtained by using our interpolation methods and those cali- 
brated by assuming one particular cosmological model (the 
ACDM model or the Riess cosmology ) with the same sample 
at z < 1 -4. This is not surprising, since the two cosmologi- 
cal models are consistent with the data of SNe la. Therefore, 
the GRB luminosity relations calibrated from the cosmologi- 
cal models should not be far from the true ones. However, we 
stress again that the luminosity relations we obtained here are 
completely cosmological model independent. 

In order to constrain the cosmological parameters, we have 
applied the calibrated relations to GRB data at high redshift. 
Since our method does not depend on a particular cosmo- 
logical model, when we calibrate the parameters of GRB 
luminosity relations, the so-called circularity problem can 
be completely avoided. We construct the GRB Hubble di- 
agram and constrain cosmological parameters by the mini- 
mum x 2 method as in SN la cosmology. From the 42 GRBs 
data (z > 1 .4) obtained by our interpolation method, we ob- 
tain n M = 0.25^;^ and n A = 0.75^$ for the flat ACDM 
model, and for the dark energy model with a constant equa- 
tion of state wo = -1.05+0 4y for a flat universe, which is con- 
sistent with the concordance model within the statistical er- 
ror. Our result suggests the the concordance model (wo = -1, 
Qm = 0.27, J7a = 0.73) which is mainly derived from obser- 
vations of SNe la at lower redshift is still consistent with the 
GRB data at higher redshift up to z = 6.6. 

For the calibration of SNe la, the luminosity relations could 
evolve with redshift, in such a way that local calibrations 
could introduce biases (Astier et al., 2006; Schaefer, 2007). 
However, Riess et al. (2007) failed to reveal the direct evi- 
dence for SN la evolution from analysis of the z > 1 sample- 
averaged spectrum. It is possible that some unknown biases 
of SN la luminosity relations may propagate into the inter- 
polation results and thus the calibration results of GRB rela- 
tions. Nevertheless, the distance moduli of SN la are obtained 
directly from observations, and therefore the observable dis- 
tance moduli of SNe la used in our interpolation method to 
calibrate the relations of GRBs are completely independent 
of any given cosmological model. 

Recently, there have been discussions of possible evolu- 
tion effects and selection bias in GRB relations. Li (2007) 
used the Amati relation as an example to test the possible 
cosmic evolution of GRBs and found that the fitting param- 
eters of the relation vary systematically and significantly with 
the mean redshift of the GRBs. However, Ghirlanda et al. 
(2008) found no sign of evolution with redshift of the Amati 
relation and the instrumental selection effects do not domi- 
nate for GRBs detected before the launch of the Swift satel- 
lite. Massaro et al. (2008) investigated the cosmological re- 
lation between the GRB energy index and the redshift, and 
presented a statistical analysis of the Amati relation search- 
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ing for possible functional biases. Tsutsui et al. (2008) re- 
ported a redshift dependence of the lag-luminosity relation in 
565 BASTE GRBs. Oguri & Takahashi (2006) and Schaefer 
(2007) discussed the gravitational lensing and Malmquist bi- 
ases of GRBs and found that the biases are small. Butler et al. 
(2007) claimed that the best-fit of the Amati relation presented 
in the Swift sample is inconsistent with the best-fit pie-Swift 
relation at > 5a significance. Butler et al. (2008) showed 
that the Ghirlanda relation is effectively independent of the 
GRB redshifts. Nevertheless, further examinations should be 
required for considering GRBs as standard candles to cosmo- 
logical use. 

Because of the small GRB sample used to constrain the 
cosmology, the present method provides no more accurate 
cosmological parameter constraints than the previous works 
where the GRB and SN la samples were used jointly. For ex- 
ample, Wang, Dai, & Zhu (2007) presented constraints on the 
cosmological parameters by combining the 69 GRB data with 
the other cosmological probes to get f2 M = 0.27 ±0.02. How- 



ever, our main point is not on improvement of the "nominal 
statistical" error; rather we emphasize that our method avoids 
the circularity problem more clearly than previous cosmol- 
ogy dependent calibration methods. Therefore, our results for 
these GRB relations are less dependent on prior cosmological 
models. 
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